Jet conversions in a quark-gluon plasma 
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Quark and gluon jets traversing through a quark-gluon plasma not only lose their energies but also 
can undergo flavor conversions. The conversion rates via the elastic q(q)g — » gq{q) and the inelastic 
qq <-> gg scatterings are evaluated in the lowest order in QCD. Including both jet energy loss and 
conversions in the expanding quark-gluon plasma produced in relativistic heavy ion collisions, we 
have found a net conversion of quark to gluon jets. This reduces the difference between the nuclear 
modification factors for quark and gluon jets in central heavy ion collisions and thus enhances the 
p/ir + and p/ir~ ratios at high transverse momentum. However, a much larger net quark to gluon jet 
conversion rate than the one given by the lowest-order QCD is needed to account for the observed 
similar ratios in central Au+Au and p+p collisions at same energy. Implications of our results are 
discussed. 
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I. INTRODUCTION 

One of the most interesting observations in central 
heavy ion collisions at the Relativistic Heavy Ion Collider 
(RHIC) is the suppressed production of hadrons with 
large transverse momentum P, This phenomenon 
has been attributed to the radiative energy loss of quark 
and gluon jets, produced from initial hard scattering of 
incoming nucleons, as they traverse through produced 
quark-gluon plasma (QGP) [1,0, [1]. Recent studies have 
shown that elastic scattering of quark and gluon jets in 
the QGP also leads to an appreciable loss of their en- 
ergies [1, 0] ■ Because of its larger color charge, a gluon 
jet is expected to lose more energy than quark and anti- 
quark jets. Since the ratio of high momentum protons 
and antiprotons to pions produced from the fragmen- 
tation of a gluon jet is much larger than that from a 
quark jet and there are more gluon than quark jets in 
proton-proton collisions, a larger gluon than quark en- 
ergy loss would lead to smaller p/ir + and p/ir~ ratios at 
high transverse momentum in central heavy ion collisions 
than in proton-proton collisions at same energy Q . This 
is in contrast to p/ir + and p/tt~ ratios at intermediate 
transverse momentum where they are enhanced in cen- 
tral heavy ion collisions as a result of quark coalescence or 
recombination [§, [l^, HH • Experimentally, data from the 
STAR collaboration have indicated, however, that p/ir + 
and p/tt~ ratios at high transverse momentum in central 
Au+Au collisions [13] approach those in p+p and d+Au 
collisions [HI]. This observation implies that the ratio of 
final quark and gluon jets at high transverse momentum 
is similar to that of initial ones. A possible mechanism 
for reducing the effect due to different quark and gluon 
jet energy losses in QGP is to allow conversions between 
quark and gluon jets via both elastic gq(q) — > q(q)g and 
inelastic qq <-*■ gg scatterings with thermal quarks and 
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gluons in the QGP. The idea of a quark jet converting 
to a gluon jet has previously been considered in Refs. flij 
for deeply inelastic scattering. In the present paper, we 
study this effect in the lowest-order in QCD. We find that 
conversions between quark and gluon jets indeed lead to 
an increase in the final number of gluon jets in central 
heavy ion collisions than the case without conversions but 
the effect is not large enough to explain the experimental 
observation. To bring the resulting p/ir + and p/-x~ ra- 
tios at high transverse momentum to the values observed 
in p+p collisions at same energy requires a factor of four 
or larger in the increase of the net quark to gluon jet 
conversion rate obtained from the lowest-order QCD. 

II. JET CONVERSION RATES IN 
QUARK-GLUON PLASMA 

The conversion rate of a quark jet to a gluon jet or vice 
versa due to two-body scattering with quarks and gluons 
in a QGP is given by the collisional width 

^ _t_ f g 2 d 3 p 2 d 3 p 3 d 3 pj 

' ' ~ 2Ei J {2ir) 3 2E 2 (2ir) 3 2E 3 (2ir) 3 2E± 
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divided by h. In the above, E\ and £3 are, respectively, 
energies of the jet before and after conversion, while Ei 
and E4 are those of quarks and gluons with degeneracy 
52 in the QGP; /(p) is the thermal distribution with + 
and — referring to gluons and quarks, respectively; and 
\M 12^34 1 2 is the squared scattering amplitude after av- 
eraging over the spins and colors of initial partons and 
summing over those of final partons and is well-known for 
both the elastic scattering q{q)g — > gq{q) and the inelas- 
tic scattering qq <-> gg. To ensure that the quark (gluon) 
jet is converted to a gluon (quark) jet in elastic scatter- 
ing, the gluon (quark) in the final state is required to 
have a larger momentum. Since both elastic and inelas- 
tic scattering are dominated by forward-peaked t— and 



2 



backward-peaked u— channel diagrams, the momentum 
of the converted jet is very close to that of the initial jet. 
Neglecting the effect due to the soft parton in a conver- 
sion process is thus justified. 
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FIG. 1: (Color online) Collisional widths for gluon to quark jet 
(upper panels) and quark to gluon jet (lower panels) conver- 
sions in QGP due to elastic scattering q(q)g — > gq(q) (dashed 
lines) and inelastic scattering qq <-> gg (dotted lines) as well 
as their sum (solid lines) as functions of temperature for dif- 
ferent quark and gluon momenta of 6 (left panels), 8 (middle 
panels), and 10 (right panels) GeV/c. 

To take into account medium effects, we include the 
thermal masses m q = m g /V3 = gT/V& [H for quarks 
and gluons in a QGP at temperature T, where g is the 
QCD coupling constant. With a s — ,g 2 /47r = 0.3, ap- 
propriate for energy scales considered here, calculated 
collisional widths for gluon to quark jet (upper panels) 
and quark to gluon jet (lower panels) conversions in a 
chemically equilibrated QGP are shown in Fig.l for jets 
of momentum 6 (left panels), 8 (middle panels), and 10 
(right panels) GeV/c. Because of larger (about a factor 
of two) quark than gluon densities in the chemically equi- 
librated QGP with thermal quark and gluon masses, con- 
tributions from elastic (dashed lines) and inelastic (dot- 
ted lines) scattering to conversion of gluon jets to quark 
jets are comparable, while inelastic scattering is more 
important than elastic scattering for quark to gluon jet 
conversion. Adding both contributions leads to a larger 
total conversion rate for the quark jet than for the gluon 
jet, particularly at high transverse momentum, as shown 
by solid lines in Fig.l. 



III. DRAG COEFFICIENTS FOR JETS IN 
QUARK-GLUON PLASMA 

As in our study of heavy quark energy loss [lq . the 
momentum degradation of a quark or gluon jet in the 
QGP depends on its drag coefficient, which is given by 
averages similar to that for the collisional width, i.e., 



where p and p' are the momenta of the jet before and 
after a collision, respectively, and the sum is over all scat- 
tering processes. Since the present study is mainly con- 
cerned with conversions between gluon and quark jets in 
the QGP, we only consider explicitly the contribution of 
two-body scattering to their energy losses and mimic the 
effect of the more important radiative energy loss by in- 
troducing a phenomenological multiplication factor Ke- 
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FIG. 2: (Color online) Drag coefficients for gluon and quark 
jets in QGP due to non-conversion (upper panels) and conver- 
sion (lower panels) two-body scattering as functions of their 
momentum in QGP at temperature T — 300 MeV (left pan- 
els) or T = 200 MeV (right panels). 

In Figl^l we show the gluon and quark drag coeffi- 
cients 7at due to non-conversion (upper panels) and 7c 
due to conversion (lower panels) scatterings as functions 
of their momentum in a QGP at temperatures T = 300 
MeV (left panels) or T = 200 MeV (right panels). It 
is seen that the drag coefficients due to non-conversion 
scattering are much larger than those due to conversion 
scattering. While the drag coefficient for the gluon jet 
is about a factor of two large than that for the quark 
jet in non-conversion scattering, they are similar in con- 
version scattering. Enhancing the calculated drag coef- 
ficients due to two-body scattering by a factor Ke = 4, 
we find that both their values and transverse momentum 
dependence are consistent with those extracted from the 
energy loss formula derived in Ref . 0] for radiative energy 
losses of quark and gluon jets in a quark-gluon plasma. 
We have also checked that with such enhanced drag co- 
efficients, the observed quenching of high transverse mo- 
mentum pions at RHIC can indeed be described. 



IV. INITIAL JET DISTRIBUTIONS AND 
HEAVY ION COLLISION DYNAMICS 

To see the effect of conversions between quark and 
gluon jets on their energy losses in QGP, we consider cen- 
tral Au+Au collisions at center-of-mass energy ^snn = 
200 GeV. The initial transverse momentum spectra of 
quark, anti-quark, and gluon jets at mid-rapidity are ob- 
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TABLE I: Parameters for initial minijet parton distributions 
given in Eq.© at midrapidity from Au+Au at -^/sjvjv = 200 
GeV. 
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670 


1.6 


7.9 
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190 


1.9 
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tained from multiplying those from PYTHIA for p+p 
collisions at same energy by the number of binary colli- 
sions (~ 960) in central Au+Au collisions. The resulting 
transverse momentum spectra of initial partonic jets can 
be parameterized by 

V(l+Pr/^r (3) 

d z p T 

with i — q,q, and g. With transverse momentum px in 
unit of GeV/c, values of the parameters Aj, Bi, and a.i 
in above equation are given in Table |H 

For the dynamics of formed QGP, we follow that of 
Ref. [ljj by assuming that it evolves boost invariantly in 
the longitudinal direction but with an accelerated trans- 
verse expansion. Specifically, its volume expands in the 
proper time r according to V(r) = ttR(t) 2 tc, where 
R(t) = Rq + ci(t — to) 2 /2 is the transverse radius with an 
initial value Ro=7 fm, ro=0.6 fm/c is the QGP forma- 
tion time, and a = 0.1c 2 /fm is the transverse accelera- 
tion. Starting with an initial temperature Tj = 350 MeV, 
time dependence of the temperature is obtained from en- 
tropy conservation, leading to the critical temperature 
T c = 175 MeV at proper time r c = 5 fm/c. 

For a quark or gluon jet moving through the QGP, the 
rate for the change of its mean transverse momentum 
(pt) is given by 

= -( 7 (p T) T)p T ) » - 7 «Pr>,T)<PT>- (4) 

In obtaining the last expression of above equation, we 
have neglected the dispersion of the jet momentum dur- 
ing its propagation and thus assumed that (p T ) « (pt) 2 - 
This is expected to be a reasonable approximation for 
high transverse momentum particles. 

Because of conversion scatterings, a quark or gluon jet 
can be converted to a gluon or quark jet with a rate 
given by corresponding collisional width. These effects 
are modeled by introducing a large number of test quark 
and gluon jets that are distributed in the transverse plane 
according to that of underlying binary nucleon-nucleon 
collisions. Their transverse momentum distribution is 
taken to be uniform with directions pointed isotropically 
in the transverse plane. A test jet of a given kind with 
certain transverse momentum is assigned a probability 
that is proportional to the corresponding jet momentum 
spectrum with the proportional constant determined by 
requiring that the sum of the probabilities for all test jets 



of this kind is equal to their total number. Motions of 
the jets are then followed via straight trajectories. Af- 
ter a small time step At, which is usually taken to be 
0.1 fm/c, a random number x is chosen and compared 
with y = TcAt/K. If x is smaller than y, then the jet is 
converted to a different kind with a reduction of its mo- 
mentum pt by 7cArpr but without changing its asso- 
ciated probability. Otherwise, the jet loses a momentum 
jnAtpt- This sequence of processes is repeated until 
the smaller of the time when the QGP phase ends and 
the time for the jet to escape the expanding QGP. The 
final momentum distributions of quark and gluon jets are 
obtained from those of test jets by multiplying with their 
associated probabilities. 

V. RESULTS 
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FIG. 3: (Color online) Nuclear modification factors for quark 
(upper lines) and gluon (lower lines) jets in central Au+Au 
collisions at ^/snn = 200 GeV as functions of momentum 
without (dash-dotted lines) or with different enhancement 
factors Kc = 1 (dotted lines), Kc = 4 (dashed lines), and 
Kc = 6 (solid lines) for conversion scattering. 

The ratio of final quark or gluon jet spectrum to its ini- 
tial spectrum, defined as its nuclear modification factor 
Raa, is shown in Figj3] Upper and lower dash-dotted 
lines are those for the quark and gluon jets using drag 
coefficients from two-body elastic scattering that are en- 
hanced by the factor Ke — ^ but without jet conversions. 
The smaller Raa for gluons than for quarks is a result of 
their larger energy loss. Although final high transverse 
momentum quark and gluon jets are initially produced 
near the surface of the QGP fireball, quark jets tend to 
originate from a region that is deeper inside the fireball 
than that for gluon jets as a result of their smaller drag 
coefficients. With the larger quark to gluon conversion 
width than that for gluon to quark conversion, it is thus 
more likely for quark jets to convert to gluon jets than for 
gluon jets to convert to quark jets. Including conversions 
between quark and gluon jets through conversion scatter- 
ings thus reduces the difference between the quark and 
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gluon Raa as shown by dotted lines in FigO 
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FIG. 4: (Color online) p/ir + (left panel) and p/ir~ (right 
panel) ratios from quark and gluon jet fragmentation in cen- 
tral Au+Au collisions at ysjvjv = 200 GeV as functions of 
momentum without (dash-dotted lines) or with different en- 
hancement factors Kc = 1 (dotted lines), Kc = 4 (dashed 
lines), and Kc = 6 (solid lines) for conversion scattering. 
Dash-dot-dotted lines correspond to p+p collisions at same 
energy. 

To obtain the proton, antiproton and pion spectra 
at high transverse momentum from those of quark and 
gluon jets, we use the AKK fragmentation function [l8l |. 
which has been shown to reproduce measured charged 
pion, proton, and antiproton spectra at high transverse 
momentum in p+p and d+Au collisions at -\/s NN = 200 
GeV. In particular, the resulting p/p ratio at high trans- 
verse momentum in p+p collisions is about 0.7 and agrees 
with the observed value. While protons and antiprotons 
are equally produced from gluon fragmentation, fragmen- 
tation of quark and antiquark jets is known to produce 
mainly protons and antiprotons, respectively p^.ll3l[l9|. 
We therefore assume that no antiprotons are produced 
from a quark jet and no p rotons are produced from an 
antiquark jet as in Ref.hJ. 

Results for p/?t + and p/-n~ ratios in central Au+Au 
collisions at ^/snn = 200 GeV are shown in the left 
and right panels Fig(H respectively. Dash-dotted lines 
are for the case without conversions between quark and 
gluon jets, and they are significantly smaller than those in 
p+p collisions shown by dash-dot-dotted lines. Including 
conversions between quark and gluon jets increases p/tt + 
and p/n~ ratios, as shown by dotted lines, but are still 
below those from p+p collisions. 

The p/tt + and p/k~ ratios are increased or the differ- 
ence between the quark and gluon Raa is reduced if the 
conversion widths shown in Fig.l are enhanced by a fac- 
tor Kc- As shown in Figs. [3] and [U with Kc = 4 (dashed 
lines), similar to the enhancement factor Ke for quark 
and gluon drag coefficients, and larger, e.g., Kq = 6 
(solid lines), the quark and gluon nuclear modification 
factors become closer and final p/n + and p/^~ ratios are 



also more similar to those in p+p collisions. 

We note that quark and gluon jet conversions have 
a larger effect on the p/ir~ ratio than the p/ir + ratio. 
This is due to the smaller number of antiquark jets than 
quark jets in heavy ion collisions and the assumption that 
quark and antiquark jets fragment only to protons and 
antiprotons, respectively. Including some admixture of 
proton and antiproton production from quark jets as well 
as from antiquark jets would reduce the difference in the 
jet conversion effect on p/ir + and p/^ ratios. 

As to the effect of using different fragmentation func- 
tions, we do not expect them to change our conclusions. 
For example, if we use instead the KKP fragmentation 
function [201 ] , which describes the pion spectrum well but 
underestimates the proton and antiproton spectra at high 
transverse momentum by about a factor of three, jet con- 
versions are still needed to obtain same p/ir + and p/7r~ 
ratios in both cental Au+Au collisions and p+p collisions 
at same energy. Of course, their values are about a factor 
of three smaller than measured ones. 



VI. DISCUSSIONS AND CONCLUSIONS 

We have studied the effect of both elastic and inelastic 
two-body scatterings of quark and gluon jets in a quark- 
gluon plasma not only on their energy loss but also on the 
conversions between them. Although inelastic two-body 
scatterings of quark- antiquark annihilation and gluon- 
gluon fusion always lead to conversions between quark 
and gluon jets, elastic two-body scatterings can be non- 
conversion scatterings if the momentum of the quark or 
gluon jet remains to be the one with a larger transverse 
momentum after the scattering. To mimic the effect of 
radiative energy loss of quark and gluon jets in the quark- 
gluon plasma, we have multiplied the drag coefficients 
calculated from two-body scattering by a factor of four. 
We have found that two-body conversion scatterings lead 
to a small net conversion of the quark jets to the gluon 
jets, resulting in slightly larger p/n + and p/tt~ ratios 
than in the case without including conversions between 
gluon and quark jets but not large enough to bring these 
ratios in central Au+Au collisions to those in p+p col- 
lisions at same energy. A large conversion enhancement 
factor of more than four, similar to that needed for the jet 
drag coefficients to describe the jet energy loss, is needed 
to explain the experimental observations. We note that 
somewhat smaller enhancement factors would be needed 
if the time for starting the interactions of initial jets with 
the QGP is earlier, such as the formation time ~ 1/pt of 
the jets, instead of 0.5 fm/c used in present study. For 
example, with a starting time of 0.05 fm/c reduces the 
enhancement factors by about 20%. 

In our study, we have assumed that the QGP pro- 
duced at RHIC is in chemical equilibrium with about 
twice many quarks and antiquarks than gluons. If the 
produced partonic matter is a pure gluon matter, result- 
ing p/n + and p/ir~ ratios in central Au+Au collisions 



5 



turn out to be slightly smaller than those from a chemi- 
cally equilibrated QGP as the gluon matter enhances the 
conversion of gluon jets to quark jets via the inelastic 
scattering gg — ► qq. On the other hand, a pure quark and 
antiquark matter gives slightly larger p/ir + and p/tt~ ra- 
tios than from a chemically equilibrated QGP as the rate 
for the inelastic conversion process qq — > gg is enhanced. 
As for a chemically equilibrated QGP, none of these two 
scenarios is able to increase the p/tt + and p/ir~ ratios 
at high transverse momentum in central Au+Au colli- 
sions to approach those in p+p collisions at same energy 
without requiring a large enhancement factor for the net 
quark to gluon jet conversion rate. 

Our results that the net quark to gluon jet conver- 
sion rate in central heavy ion collisions at RHIC is much 
larger than that given by the lowest-order in QCD may 
not be surprising as the jet drag coefficients from the 
lowest-order QCD need to be enhanced by a similar fac- 
tor to describe their energy losses in QGP. Also, previous 
studies using the multi-phase transport (AMPT) model 
[2lj . that includes only two-body scattering among par- 
tons, have shown that a much larger parton scattering 
cross section than that given by the lowest-order QCD is 
needed to describe many other experimental observations 
at RHIC such as the the larg e elliptic flows of hadrons 
made of light quarks [13, l23f or heavy quarks [24[ and 
the two-pion correlation functions [251 ] - The large en- 
hancement factor over the lowest-order QCD results can 



be considered as an effective parameter for taking into 
account effects not considered in present study such as 
higher-order contributions and multi-body scattering. A 
preliminary study of jet conversions via higher-order ra- 
diative scatterings with thermal partons shows that they 
lead to similar quark and gluon jet conversion widths 
as the lowest-order two-body scatterings considered in 
present study. Also, using the QCD coupling constant 
at finite temperature as extracted from the lattice QCD 
studies [26| . which is about a factor of two larger than 
that from the perturbative QCD, further increases the 
quark and gluon jet conversion widths, resulting in a con- 
version enhancement factor close to the needed one. Our 
results thus could be another indication for the strongly 
coupled QGP that has been produced in central heavy 
ion collisions at RHIC [13]. 
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